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In this paper we investigate x-ray and 7-ray propagation in crystals having a con- 

OO 

stant strain gradient and flat or cylindrical surfaces. When a displacement field is 
present, we solve the Takagi-Taupin equations either by the Riemann- Green method 
or by a numerical method. We apply the results to study the operation of a double- 
crystal Laue-Laue diffractometer having a flat collimating crystal followed by a bent 
analyzer crystal. In particular, we examine the effect of the analyzer strain on the 
location of the diffraction peaks in the dispersive and non-dispersive setup, thus con- 
firming our previously reported peak-location as being set only by the diffracting plane 
spacing on the analyzer entrance surface. 
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1. Introduction 

The main purpose of this work is the to verify the soundness of what we had previously 
published (Mana et al, 2004). On that occasion, we had studied how a constant 
strain gradient in the rotating crystal of a double-crystal diffractometer affects the 
instrument operation. The result of that investigation was that the position of the 
Bragg peaks depends only on the diffracting-plane spacing on the crystal entrance- 
surface. However, to simplify the problem, we had then assumed all the crystal surfaces 
rigorously flat. Since subsequent experiments, using cylindrically bent crystals and 
both x- and 7-rays to test that rather surprising result, delivered contradictory results 
still under examination (Kessler 2007, Massa et al., 2005), we were urged to investigate 
propagation in bent crystals in more detail. An additional reason was the interest in 
efficient Laue-Laue bent crystal diffractometers for 7-ray spectroscopy of nuclei having 
a very high thermal neutron cross-section (Materna et al., 2006). 

For these reasons, we extend here our previous analysis by taking account of the cur- 
vature of the crystal surfaces and by simulating diffraction in both the non-dispersive 
and dispersive setups. In section 2 we solve, by the Riemann-Green method in carte- 
sian coordinates, the Takagi-Taupin equations for the propagation of x- and 7-rays 
in bent crystals; we indicate also how the crystal surfaces are modelled, what choice 
of the reference perfect-crystal we have adopted and what kinds of distortion of the 
cylindrically bent crystal we have considered. Sections 3.1 and 3.2 deal with two cases 
when the crystal surfaces, on which the initial conditions have to be assigned, are flat 
or cylindrical, respectively. While in the first case the solutions are known (Authier 
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& Simon 1968, Mana & Palmisano 2004), and are re-examined here to illustrate our 
formalism, we are not aware of solutions when the crystal surface is a cylinder. In 
fact, in the literature, the case of a curved crystal surface in the macroscopic sense 
is just hinted in Takagi (1969) by means of a curvilinear coordinates. Subsequently, 
Olekhnovitch & Olekhnovitch (1980) carried out the calculation of the profile function 
of the scattering curve for a crystal in the form of a cylinder the size of which does 
not exceed the extinction length. Later, Thorkildsen & Larsen (1998) observed that 
it is in principle possible to obtain analytical expressions for the primary extinction 
factor in perfect crystals having a circular diffraction plane. In section 4 we apply 
our results to a double diffractometer and, to corroborate them, in section 5 we solve 
numerically the Takagi- Taupin equations in polar coordinates for different geometrical 
and physical parameters. 
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2. Takagi-Taupin equations for distorted crystals 

In order to study x-ray propagation through a distorted crystal, we shall apply the 
Takagi-Taupin equations (Takagi 1962, Takagi 1969, Taupin 1964, Autier 2001, Mana 
& Montanari 2004) in the two-wave approximation of the dynamical theory of x-ray 
diffraction. Since only cylindrical geometries will be considered, we shall use a two- 
dimension model with a reference frame having the x and z axes lying in the reflection 
plane. Lattice distortion is described by the displacement field u(x, z), which gives the 
difference between the actual distorted lattice and a virtual perfect-lattice identified 
by the reciprocal vector ho, which will be chosen according to our convenience. Hence, 
we anchor the reference frame to it and set the x axis parallel or antiparallel to ho. 
By using directional derivatives, the Takagi-Taupin equations can be written as 



= ~^ Do + — Dh (la) 

= ^D h+ ^D 0+ d -^D h . (lb) 
ash 2 2 ash 

In equations (la) and (lb), D Q and Z\ are slowly varying amplitudes of the Ewald 
expansion 



D = D exp (iK • r) + D h exp [i(K h ■ r - h • u)] (2) 
of the dielectric displacement vector D = Dy for the a polarization, 

s = K Q = — xsin + z cos 6*b (3a) 
§h = Kh = x sin + z cos ^b (3b) 

are the unit propagation vectors, ho = x2ifsin6>B, Kh = K Q + ho, K = ||K || = 
||Kh|| = 27r^/c is the modulus of the wave number vector of the incoming radiation 
(with frequency v), 6b is the Bragg angle (with a sign), and the complex parameters Xo, 
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Xh, an d X-h are the Fourier components of electric susceptibility; in our case x-h = Xh- 
We shall consider an infinite crystal slab the surfaces of which are smooth curves, 
r : [n, r 2 ] -► R 2 and S : [cti, a 2 ] -► R 2 , with ||dr/dr|| and ||dS/dcr || e]0, l/tan6» B [, 
where the last constraint is necessary to have a well-posed problem. By imposing that 
at each point of the entrance surface the o-component of the Ewald expansion (2) 
is equal to the incoming wave D c (x,z) = A(x, z) exp(iK c r), the initial conditions to 
solve equations (la, lb) are therefore 

D (x,z)\ r = $(s,z)| r (4a) 
D h {x,z)\ r = 0, (4b) 

where 3>(x, z) = A(x, z) exp [i(K c — K )r] . 

Equations (la, lb) can be simplified by introducing two new amplitudes, d Q and dh, 
defined as 

/ .Kxo s + s h \ 
do,H = exp^-i— rT§ - 7§h -rjD , h . (5) 

By substituting equations (5) for D Q ^ into equations (la, lb), we obtain 



— l 



.ddo 
ds Q 
. dd h 
ds h 



2 



dh 



g(ho-u) 
2 as h 



with the initial conditions 



d (x,z)\ r = exp ( — - ^° j- ^ • r ) 3>(x, z) 



2 1 + s • s h 



dh(x,z) 



= . 



(6a) 
(6b) 



(7a) 
(7b) 



As in the cylindrical crystals the geometry is somewhat elaborate, in the following 
subsections we illustrate the main underlying assumptions. 
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2.1. Crystal surfaces 

Let us now exemplify how cylindrical surfaces are modelled. The curvature centre 
can be located either on the source side or on the opposite, with Rq and Rq + T 
(T denoting the crystal thickness) being the curvature radii of the entrance surface, 
respectively. A crystals having cylindrical surfaces will be called concave when its 
concavity is towards the source, convex otherwise. 

In the concave case, by locating the reference-frame origin at the curvature centre, 
the parametric components of the surfaces are 

f r x (r) = i? sinr . . 

\ T z (t)=R cost { ' 

and 

f £ x (<7) = (Rq + T) sin ct . . 

\ E*(a) = (R + T) cos a ' { ' 

where r G [t± , T2] and cr G [ui , C2] . The inward and outward normals to the entrance 

and exit surfaces, respectively, are 

I (n r ) (r)=sinr 
( (nr)z(r) = cost 

and 

f (n s ) x (*) = sma 
{ (n s ) 2 (cr) = coscj 

In the convex case, by locating the reference-frame origin again at the curvature centre, 
the surface components and their inward and outward normals are 

/ r x (r) = (Ro + T)sinr . . 

\ T z (t) = (R + T)co S t ' 

J Z x (cr) = R s'ma . . 

I S 2 ((7) =^0 003(7 ' 1 j 



IUCr macros version 2.1/31: 2007/05/15 



/ (ftr)*(r) = 
1 (ftr),(r) = 



— sinr 

— cos r ' 



(14) 



and 

j (n s ) x (a) = -sino- 



- • (15) 



Crystal rotation 

With a rotated crystal, since the reference perfect-lattice can be chosen at our con- 
venience, we make the virtual lattice immovable, no matter what the crystal rotation 
might be; consequently, a rotation is nothing that a very special kind of distortion. 
Hence, our Ansatze to study x-ray propagation in a rotated crystal are: firstly, ho is 
independent of rotations, secondly, the first-order rotation term 

u rot (x, z) = [zsina + x(cos a — l)]x ps azx, (16) 

where a is the rotation angle (clockwise oriented), must be added to the crystal dis- 
tortion and, thirdly, the crystal surfaces must be mapped into the new lines (1^,1^) = 
(T x cos a+T z sin a, — T x sin a+T z cos a) and (T,' x , T,' z ) = (S x cos a+S z sin a, — T, x sin a+ 
S^cosa). As long as a < 1, we shall consider cos a « 1 and sin a w 0, so that the 
crystal surfaces are left unchanged by rotations. 

2.3. Lattice distortion 

We shall consider two distortions describing cylindrically bent crystals; both, fan- 
down, 

(1)/ \ _ x ( z ~ z m) 



and fan-up, 



«4 lf (*,*)= „ +r\ (17a) 



v x 2 \x,z) = ^-^, (17b) 



are pure displacements in the x direction. We have chosen the x axis origin and the 
reference perfect-lattice so that, for the non-rotated crystal, u(x = 0, z m ) = and 
ho = h(x,z = z m ), where h is the local reciprocal vector of the distorted lattice. 
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In equation (17a), the displacement is such that all the lattice planes are directed 
towards a point at distance Ro from the entrance surface and the points at z = z m 
are undisplaced (left upper part of Fig. 1). In equation (17b), the lattice planes are 
directed towards a point at distance Rq from the exit surface. For a rotated crystal, 
the complete displacement field, including both equations (17a, 17b) and equation 
(16), is 

uf{x,z) = v${x,z) + u^{x,z) 

= g^^^ + az, (18) 

where g(l) = 1, g{2) = —1 and the approximation v x % \x, z) cos a « Vx\x,z) is 
applied. 



3. Propagation in distorted crystals 

Equations (17a) and (17b) approximate the displacement field, characterized by a 
constant strain gradient, in crystals having their surfaces fiat or cylindrical. 



3.1. Flat crystal surfaces 

In the simplest case of fiat external crystal surfaces, the Takagi-Taupin equations 
are 



sin #b 



+ cos6> B 
ox oz 



<9do } . Kx-h ,(i) 



*inf) , rn *f> 9d h _- K Xh,(i) 



+ \2K sin <9 B 
with the initial conditions 



9{i) p , Zm sin B + (g(i) — — V a ) cos 6 B 

-Ko + z m \ -Ko + z m 



(19a) 



4° (19b) 



d«(x,0) = t>(x) 
d£\x,0) = 0, 



(20a) 
(20b) 
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where we assumed that the external crystal surfaces are orthogonal to s Q + §h (sym- 
metrical Laue geometry) and we located z = on the entrance surface. 

The coupled equations (19a, 19b) can be simplified by a change of dependent vari- 

~ (i) ~ (i) 

ables. Let us introduce the two unknowns Do and defined by the expression 

^l 1 = exp(-i/W2^sin^ B )< ) h , (21) 
where the function f( l \x,z) is defined as 

f (i) (x,z) = g(i)r—^- 7 + - V) sin ^BCos0 B - zj- (^- r + — 

Ro + Zm 4 Vsmt'e cos6'b/ 2 \sm0B cos 0b 

1 / x z \ 

+ a- -— - + — cos 6» B . 

2 \sin^B cospb/ 

If we observe that 

d 



(^B^+COS0 B ^)/ ( *W) 



sin# B g(i) { * * m) +cosg B (//;'),, • n) (2:3) 

-Ko + z m \ Kq + z m 



and 



(" Sin ° B Tx + C ° S ° B §z) /W (X ' Z) = ° ' (24) 
substitution of equation (21) into equations (19a, 19b) gives the Takagi-Taupin equa- 
tions in the unperturbed form 

sin0B ^ + cos ^ B<9 lr = iK ?^ l) ■ (25b) 

As a consequence of equation (21) the new initial conditions on the entrance surface 
z = are 

D®(x,0) = ex.p[-if (i \x,0)2Ksm9 B }^(x) (26a) 
I)£\x,0) = 0. (26b) 
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By the Riemann- Green method (Authier & Simon, 1968, Takagi, 1969, Sommerfeld, 
1964, Palmisano & Zosi, 2005) we can find a solution by quadrature for the system 
(25a, 25b) with initial conditions (26a, 26b) 

D%(x,z) = / +O °G Oih (x-x , ,z)exp[-i/W(x , ,O)2^sin0B]^(^)d^/ (27) 

J — oo 

where the kernels G Q and are 

G (x, z) = <5(x + ztan^e) 
K 

~ a\ ■ a i \/XhX-h H(z|tanfl B | + x)H(z| tan fl B | - x) 
4 sm tf B 



' z tan 6> B — x / K 



z tan 9b + x \2 \ sin # B 



Ji( 2|s .^ B| V / XhX^rv / - z2tail2g B-^) (28) 



and 



Gh(x,z) = ^^H(z| tan 6* B | + x)H(z| tan # B | — x) 

4 | sm 9b I 

x Jo ( 2 \snl9B\ ^^^ z2tan2dB ~ x2 ) ' (29) 
In equations (28) and (29) H(z) is the Heaviside function and Jo(z) and Ji(z) are the 
Bessel functions of the first kind and order 0, 1, respectively. Therefore the solution 
to Takagi- Taupin equations (19a, 19b) with initial conditions (20a, 20b) is 

d%{x,z) = exp(i/«(x,z)2Ksin# B ) 

x / G Q:h (x- x',z)exp[-if< l \x',0)2Ksin9B]^(x')dx' . (30) 

We see from equation (30) that the effect of a constant strain gradient has been reduced 
to a similarity transformation of the G Q and Gh kernels (Mana & Palmisano 2004). 
Equations (5) and (30) show that the intensities of the transmitted and diffracted 
beams Do\x, z) and D^\x, z) on the exit surface z = T are 

f + OO 



= / D^ h (x,T) cos0 B dx 

J — oo 

/+oo 
ex V [-KQ( X o)T/cos9B} 
-oo 



' — oo 
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/+oo 
G Q , h (x - x', T) exp[-i/ w (x, Q)2Ksin9 B ]^{x) dx 
-oo 



cos #b dx , 



(31) 



where 9(xo) is the imaginary part of % . Eventually, substitution of equation (22) into 
equation (31) gives 



+oo 



exp[-K9 ; (xo)r/ cos 6> B ] 



/+oo 
G , h (x-x',T) 
-oo 



x exp | -i2i-C sin B 
x cos #b dx . 



cos Ob 
Rq + z„ 



1 ^ 



12 



sin#B 2 V cos#b sin#B 



a 



<f>(x')dx' 
(32) 



Equation (32) gives the rocking curves /^(a) when the crystal is distorted by the 
displacement field (17a) or (17b), the external crystal surfaces are fiat and <£(x) is the 
complex field amplitude of a generic incoming beam. Additionally, and generally, the 
effect of the displacement fields (17a, 17b) or (18) on the intensity in equation (32) 
is seen to consist of a phase-redefinition of the initial condition 3>(x). Evaluation of 
equation (32) in the limit with Rq tending to +oo gives the rocking curve of a perfect 
analyzer crystal. With the variable change a = + a' , where is 



a 



(0 - 



9{i) 



— — tan B 

(x, 0) tan Ob , 



dx 



equation (32) can be reduced to the simpler form 



(33a) 
(33b) 



/«(««+«') 



+oo 



exp[-KQ(xo)T/ cos B ] 



/+oo 
G , h (x-x',T) 
-oo 



x exp < -\2K sin B 



cos #b 

Rq + Z n 



. 1 / Ro + z m i 

+ -x — ; — - — a 



- O — 

4 9{l) sm0 B ' T sin^B 



3>(x')dx' 



cos Ob dx . 



(34) 
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3.2. Cylindrical surfaces 

By application of the Riemann- Green method, the solutions of system (19a, 19b), 
with the initial conditions (7a) and (7b), are the flux integral 



exp[ifW(x,z)2Ksm9 B } 

J G , h (x -T x ,z- T z ) exp[-i/«(r)2^sin 9 B ] 
. Kxo s + s h 



x exp 



(T x x,T z z) 



COS 6>B 



(35) 



2 1 + s • s h 

which generalizes equation (30). The unit vector fir is the inward normal to the 
r : [ti,T2] — > R 2 surface, T x and T z are the surface parametric components, dr is a 
shorthand form for ||dr/dr|| dr, and < ||dr/dr|| < l/tanf9 B Vr G [t 1 ,t 2 \. Let n s 
be the outward normal to the exit surface £ . From equations (5), (35) and (22) the 
intensity of the forward transmitted and diffracted beams are 



= jf exp (-if9f(xo) -J^- ) ^ Go.hCEa, - r x , S 2 - T 2 ) expj - i2K sin B 



sin 6>b cos #b 

#0 + 2m 



+ 



r, 



4 Vsin^B cos#b 



1 



r. 



+ 



r 2 



X exp ( — i — - — 7777T" ) $(1^,1^) — dr 



2 \sin#B cos^b 



s h • n s dE . 



cos #b sin #b 



c 



2 COS 6*B / cos #B 

In the concave case, the equation (36) with use of the a^ l \ and the definitions (33a) 
and (8-11), can be rewritten as 

(i?o + T) cos a — Rq~ 



cos 8b 



jg(a< 1 > + o') = Texp 

/ G 0i h((Ro + T) sin a — Rq sin r, (i?o + ^) cos <* — Ro cos r) 

. „ sin^BCOs^B^ / sinr cost— 1 
-i2ifsin6> B „ — -Ro ( -^r- + 



Ro + z n 



sin #b cos 9b 



Rq / sinr cosr-l \ 1 i? + z m / 
51 j 4 Un0B cos^b / 2 sin# B 



x exp — i 



Rq 



cos r — 1 
cos 9b 



$(i?osinr, Rq(cost — 1)) 
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R dr 



cos(6* B - a)(R + T)da. 
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(37) 



cos 6b 

In the convex case, the equation (36), by following the same procedure, that is, with 
the use of definitions (33a) and (12-15), can be rewritten as 

Rq + T — Rq cos a ' 



cos #b 



^ 2 h(« (2) + «')= Texp -K$S( Xo )- 

/ G Q> h{Ro sin a — (Rq + T) sinr, Rq cos a — (Rq + T) 

Jti 



COS T 



x exp < —\2K sin #b cos 0b 



Rq + T ( sin r cos r — 1 



+ 



-Ro + -2m \sin 9b cos 0b 



x 



9(2)- 



Rq + T f sin r cos r — 1 \ 1 Rq + z m / 



+ 



+ 



x exp 

cos(6/r — t 
x 1 

cos 9b 



sin#B cos^b / 2 sin#B 
cos r — 1 



Q: 



■^(Ro + T) 

2 cos pb 



$((i?o + T) sinr, (i? + T)(l - cost)) 



(i?o + r)dr 



cos(0b + o-)i?odcr . 



(38) 



4. Double crystal diffractometer 

Figs. 1 and 2 show the Laue-Laue diffractometer in both the non-dispersive and disper- 
sive setups. We consider a fiat collimating perfect-crystal and a monochromatic point 
source located in (Lsin#B, — Lcos#b), where L is the distance between the source and 
the entrance point of the collimating crystal. We have chosen the vector ho of the 
analyzer reference perfect-lattice equal to the vector ho collimating crystal; therefore, 
these two lattices have the same spacing and, when a = 0, they are parallel. If the 
two crystals have the same lattice spacing, the analyzer rotation a between the non- 
dispersive and dispersive geometries is equal to 29b] a different rotation is related to 
a different lattice spacing in the collimating and analyzer crystals. 
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4-1. Rocking curves 

Let us confine our study to the reflected beam. When we examine the equation (34) , 
we see that, if 

$(ar) = , (39) 

where <J>(x) is the amplitude of the o-component of the external field on the entrance 
surface of the analyzer, the intensity profile of I^\a^ + a') is invariant under the 
a' T^- —a' exchange, as can be easily checked by the double substitution x' = —x' , x = 
—x . This means that I^\a) has a vertical symmetry-axis passing through aW. Let us 
note that equation (39) is fulfilled if the amplitude of the o-component of the incoming 
beam is an even function. 

The formula o^/ tan 0b = —du$/dx, obtained from equation (33b), is our formu- 
lation of the AO/ tan6*B = —Ad/d Bragg's law. It is to be noted that, if the entrance- 
surface displacement field is zero, i.e., when z m = 0, then aW = and there is no shift 
of the reflection peak with respect to the perfect-crystal case. 

With crystals having cylindrical surfaces, as long as the approximations cost 
cos(#b + t) ~ cos(#b — r) « cos(#b — f ) ~ cos(0b + a) « 1 in equation (37) and in 
equation (38) are valid, and the symmetry requirement in equation (39) is satisfied, 
the same conclusion holds. In fact I^\a), now given by equation (37) or equation 
(38) and depending on the concave or convex case, has a symmetry axis passing 
through = g(i) tan #b z m / (Ro + z m ). The approximations mean that we have small 
Bragg angles, great curvature radii and the profile of the wave from the collimating 
crystal is not altered by the geometry of the analyzer crystal. In Mana et al. (2004), 
it is reported that, in the presence of a constant strain gradient in the analyzer, 
the Laue-Laue rocking curve is shifted by (Ad/d) tan Ob, where the lattice strain is 
evaluated on the crystal surface. However, this peak shift is not easily measurable. 
Since the analyzer rotation between dispersive and non-dispersive reflection peaks 
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is an experimentally observable quantity, we give now the relevant equation. In our 
formalism, the o-component of the analyzer crystal field co-propagates or counter- 
propagates with respect to x axis according to the 9b sign; therefore, the exchange 
between the dispersive and non-dispersive geometries corresponds to the substitution 
of —9b for 9b- If we observe that the non-dispersive rocking curve peaks when a = 
9b), where 9b is the Bragg angle for the collimating crystal (Fig. 1), and we 
observe as well that the dispersive one peaks when a = 29b + sW(9b) (Fig. 2), the 
sought formula is 



29 B + a^(e B )-aW(-9 B ) 



(0/ 



9 B + a^{9 B ) 



9b 



dx 



(x, 0) tan#B 



(40a) 
(40b) 
(40c) 



where (du^ /dx)(x, 0) is the lattice strain on the analyzer entrance surface. 



5. Numerical simulation 

In order to validate our previous results we have also studied the diffractometer oper- 
ation by solving the Takagi-Taupin equations numerically. Again we consider the 
two distinct cases of fiat or cylindrical analyzer surfaces. The collimating crystal is 
a parallel-sided silicon slab limited by two surfaces orthogonal to the (220) Bragg 
planes. The x- or 7-ray source illuminates the collimating crystal by a monochromatic 

cylindrical wave D e (x, z) = *&(x, z) exp(zK • r), where 

[x-(w/2)] 8 [x + (w/2)Y 



*(x,0) = < 



A(xY- 




(w/2) 



16 



if - w/2 < x < w/2 
otherwise 



A(s) 



1 



exp 1- 



, 2-7T COS 2 #B 2 



-X 



(41) 



(42) 



4ttL "' r ^ A 2L 

A = hc/E (with E the photon energy) is the wavelength, w/2 is the half- width on the 
entrance slit, and the Bragg angle is positive. In the numerical simulation, we consid- 



IUCr macros version 2.1/31: 2007/05/15 



16 

ered the two sets of parameters shown in Table 1; the lower energy value refers to the 
experimental setup described in Mana et al. (2004), the upper in Massa et al. (2006). 
We have considered silicon (220) Bragg planes and have taken the values of the dielec- 
tric susceptibilities from the Sergey Stepanov's X-Ray Server http:/ / sergey.gmca.aps.anl.gov . 

As a first step we solved equations (19a, 19b) for a perfect crystal with boundary con- 
ditions (41). Subsequently, the reflected beam, which we shall indicate by D^' co \x, T), 
was free-propagated rigidly from the collimating crystal onto the entrance surface of 
the analyzer crystal. 

5.1. Flat analyzer surfaces 

In this case, the propagation equations of the analyzer fields Do^' ana '(x, z) and 
D^' ana (x, z) are (19a, 19b), where u£\x,z) is the displacement field in equations 
(17a) and (17b), the Bragg angle is negative, and the initial field values are 

D«' ana (x,0) = D^' co \x,T) (43a) 
Df ),ana (x,0) = 0. (43b) 

We have calculated the numerical solutions with the aid of MATHEMATICA (Versions 
5.2 and 6.1, Wolfram Res. Inc.); we have obtained the same results (to within 0.1%) by 
performing the integration in equation (30) numerically. In Figs. 3 and 4 we show the 
rocking curves when the analyzer distortion corresponds to the fan-down and fan-up 
cases, in both the non-dispersive and dispersive geometries; the peak-shifts agree with 
the values predicted by equation (33a). 

Two cases, calculated numerically and according to equation (33a), are compared 
in Fig. 5, showing that the maximum difference in Aa amounts to a few parts per 
10~ 7 . Figs. 3 and 4 require a few comments. Firstly, in contrast to experimental 
observations, there is not spreading in Fig. 4; the reason is that we assumed the 
source monochromatic and we did not integrate over its linewidth. 
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In a non-dispersive geometry, with a perfect-crystal analyzer having the same lattice 
spacing as the collimating crystal, all rocking curves peak when a = 0, no matter 
which the wavelength might be. On the contrary, in a dispersive geometry, the rocking 
curves peak when a = 29 B {\) thus giving rise to a convolution integral. Secondly, we 
did not give particular attention to the crystal field intensities; therefore, the relative 
intensities of the 17 keV and 184 keV plots are meaningless and the two curves are 
not comparable. Finally, the dotted lines show only the central part of the extremely 
wide 17 keV curve. 

5.2. Cylindrical analyzer surfaces 

In this case, we rewrote the Takagi-Taupin equations (19a, 19b) in polar coordinates, 
p = V x 2 + z 2 and ip = arctan(z/x), the reference-frame origin being in the centre of 
curvature of the entrance surface. Hence, 



... , ^/Jfi- sin(g B + ^) rf'^ KX-h nW ,ana 

C os{e B + v)— d^ = '— D * ( 44a ) 



cos { e B - + -^-^t' ana - 

op p dip 2 
+i2Ksm9 B I cos (0 B - <p)-qJ- + " \ D h ( 44b ) 

the x-component of the displacement field, including the rotation term, is 

(i) , , . (pcostp - R - z m \ 

"idlftW = psmip[ — + g(i)apcosp (45) 

V Kq + z m J 

and the boundary conditions for the concave case are 

£>W' ana (i?o,T) = D { ^ co \R G smT + R G {1 - cos t) tan fl B , T) (46a) 

£>« ,ana (i2o,r) = 0, (46b) 

and 

D«' ana (i?o + T,t) = D^' CO \(R + T)sinr - (R + T)(l - cos r) tan 6> B , T4 7a ) 
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D^' ana (R + T,T) = 0, (47b) 

for the convex case. Here, too, the numerical solution of equations (44a) and (44b) 
agrees with the numerical integration of equation (35). 

Figs. 6 and 7 show the rocking curves for a concave and for a convex crystal, in the 
non-dispersive and dispersive geometries. Also in this case the peak shifts agree with 
the values obtained from equation (33a). By comparing the profiles in Figs. 3 and 6 
we observe a slight effect depending on the surface curvature which, anyway, does not 
alter their symmetry; furthermore, the wings of the upper curves do not appear in the 
range shown; the same can be concluded from Figs. 4 and 7. 

Additionally, the figures exemplify that peak shifts and, consequently, Aa, are inde- 
pendent of the surface geometry (flat or cylindrical, concave or convex), but they 
depend only on the entrance-surface strain Ad/d. This confirms that the flat-surface 
approximation used in Mana et al. (2004) was admissible. To check further our numer- 
ical computations, we examined also two auxiliary cases, when the Bragg planes are 
simply either contracted or expanded and when the Bragg planes are not distorted. 
We have also carried out many numerical simulations with different entrance-slit aper- 
tures; the above conclusion were always confirmed. 

6. Conclusions 

We have studied x- and 7-ray propagation in flat and cylindrically bent crystals. We 
have used the relevant results to describe the operation of a Laue-Laue diffractome- 
ter consisting of a flat collimating crystal and a bent analyzer crystal and we have 
extended the results of a previous investigation of ours (Mana et al., 2004). We have 
described the distortion charaterized by a constant strain gradient in crystals having 
flat or cylindrical surfaces. 

In both cases, in addition to numerical simulations, we have also given exact solu- 
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tions of the Takagi-Taupin equations in the form of Riemann-Green integrals. We have 
confirmed both analytical and numerical results that the rocking curve shift does not 
depend on the shape of the analyzer surface, but only on the lattice strain on the 
entrance surface. 

Since the validity of these solutions - via convolution integrals - is not limited to 
fiat and cylindrical surfaces, we can extend such solutions either to the case when 
also the collimating crystal is cylindrically bent, or to the case when the effect of the 
surface roughness in x- and 7-ray diffractometry and interferometry is not negligible. 
Additionally, numerical simulations open the way to a better understanding of the 
operation of bent-crystal diffractometers, in particular when finite element solutions of 
the elasticity equations are integrated into the Takagi-Taupin equations thus allowing 
us to characterize the relevant lattice strains of the diffractometer crystals. 

This work was supported by the Museo Storico della Fisica e Centro Studi e Ricerche 
"Enrico Fermi", Rome, the Regione Piemonte and the Compagnia San Paolo, Turin. 
The Institut Laue-Langevin is gratefully acknowledged for beam time allocation and 
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Table 1. Parameters used for different energies 
ref. T (mm) L (m) w/2 (/im) E (keV) R a (m) 

Mana et al. (2004) 0.5 1 50 17 160 

Massa ei al. (2006) 2.5 16 500 184 697 




Fig. 1. Layout of a two-crystal diffractometer in non-dispersive geometry. On the left, 
the analyzer crystal surface is flat and on the right it is cylindrically bent. The 
symbol I is the incident beam, T and R are the beams transmitted and reflected 
by the collimating crystal, RT and RR are the beams transmitted and reflected by 
the analyzer crystal. The unit vectors s h are defined in equations (3a) and (3b). 
The angle of incidence 9 is positive {6 = 6-q) on the collimating crystal and negative 
(9 = —6b) on the analyzer crystal. 
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Fig. 2. Layout of a two-crystal diffractometer in non-dispersive and dispersive geome- 
tries. The symbol I is the incident beam, R is the beam reflected by the collimating 
crystal, RR nd and RR d are the double reflected outgoing beams, the superscripts 
"nd" and "d" refer to the non-dispersive and dispersive setup. The unit vectors s 0j h 
are defined in equations (3a) and (3b). The angle of incidence 9 on the collimat- 
ing crystal is positive {9 = 9b)', the angle of incidence 9 on the analyzer crystal is 
negative for the non-dispersive setup and positive for the dispersive setup. 
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Fig. 3. Non-dispersive (left) and dispersive (right) rocking curves for a fiat Si analyzer 
with a fan-down distortion. Upper and lower curves refer to the energies of 17 keV 
and 184 keV, respectively. Bragg planes are (220), the relevant parameter values 
are given in Table 1. 
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Fig. 4. Non-dispersive (left) and dispersive (right) rocking curves for a fiat Si analyzer 
with a fan-up distortion. Upper and lower curves refer to the energies of 17 keV 
and 184 keV, respectively. Bragg planes are (220), the relevant parameter values 
are given in Table 1. 
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Fig. 5. Analyzer rotation Aa from non-dispersive to dispersive geometry as a function 
of curvature 1/R for T = 1.4 mm (t), T = 2.5 mm (★) when E = 184 keV and 
z m = T/2. Positive and negative values of 1/R refer to the concave and convex 
cases, respectively, and have been calculated numerically. 
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Fig. 6. Non-dispersive (left) and dispersive (right) rocking curves for a concave Si 
analyzer with a fan-down distortion. Upper and lower curves refer to the energies 
of 17 keV and 184 keV, respectively. Bragg planes are (220), the relevant parameter 
values are given in Table 1. 
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Fig. 7. Non-dispersive (left) and dispersive (right) rocking curves for a convex Si 
analyzer with a fan-down distortion. Upper and lower curves refer to the energies 
of 17 keV and 184 keV, respectively. Bragg planes are (220), the relevant parameter 
values are given in Table 1. 
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